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^«0 , Abstract We calculate the low temperature corrections to the free energy for a 

sphere in front of a plane. First, the scalar field obeying Dirichet or Neumann bound- 
ary conditions is considered. Second, the electromagnetic field is studied, the sphere 

r^ ' being perfectly conducting and being a dielectric ball with both, constant permittiv- 

ity and permittivity of the plasma model. 



1 Introduction 



During the past decade significant attention was payed to the Casimir effect at finite 
temperature. The interest was triggered by the desire to measure the temperature 
dependent part of the force (see IIKMM09I , section 4.D, for a review) and by the 
conceptual problems arising for T — )- with some thermodynamic quantities (see 
IIKMM09I . section 2.D). Recent interest came also from the interplay between tem- 
perature and geometry investigated in IGWIOI IWGIOI using world line methods. 
For a scalar field with Dirichlet boundary conditions on the interacting surfaces, a 

^^ ■ generic behavior ~ T'^ of the temperature dependent part of the force was found and 

attributed to open geometry. 

In our previous paper IIBP09I we investigated the interaction between a ball and 
a plane for both, the scalar and the electromagnetic fields. We used the exact func- 

^^ ■ tional determinant method (also called scattering approach or 'TGTG' -formula) and 

H I focused on the limit of small separation. We showed that the Proximity Force Ap- 

proximation (PFA) is reproduced exactly for medium and high temperature. For low 
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temperature, the leading order of the free energy is the vacuum energy (for which 
the PFA is reproduced, of course) and the temperature dependent part is a small 
addendum for which the PFA does not hold. 

In the present paper we calculate the leading behavior of the temperature depen- 
dent part Aj^ of the free energy and that of the force, Ajf, for T ^ 0. Again, we 
consider the scalar and the electromagnetic fields. For the latter we allow for disper- 
sion including fixed permittivity and permittivity following from the plasma model 
(see below, eq. (|30]l). We use the functional determinant method and truncate the 
orbital momentum sum at some finite /„,. It turns out that the limit of /,„ -^ °° shows 
a sensitive dependence on the separation a for a ^ 0. Due to the space limitations 
we will be quite brief concerning technical details. These will be given in a separate 
paper. 

We mention that the interplay of geometry and temperature in sphere-plane ge- 
ometry was also studied in ICDMNLRIOII and related papers using the scattering 
approach. In a numerical analysis the ambient temperature was fixed, while the ra- 
dius of the sphere and plasma frequency varied. 

Throughout the paper we use units with h = c = 1 . 



2 The free energy at finite temperature 

In this section we collect the basic formulas for the free energy in functional deter- 
minant representation at finite temperature. We follow closely the notations used in 
IIBP09I ■ At finite temperature, the free energy is given by 

^ = 7 I Trln(l-M(<^„)), (1) 

where £,„ = InnT are the Matsubara frequencies. The matrix M results from the 
scattering on the sphere and will be described below together with the meaning of 
the trace. The sum over the Matsubara frequencies can be transformed into integrals 
using the well known Abel-Plana formula. The free energy separates into two pieces, 

^ = £o + 4r^, (2) 

where 

£o = ^/_"§Trln(l-M(^)) (3) 

is the vacuum energy, i.e., the free energy at zero temperature, and 

4r=^-^/_"^«r(<^)/Tr[ln(l-M(/^))-ln(l-M(-/^))] (4) 

is the temperature dependent part of the free energy containing the Boltzmann factor 

nr(^)=l/(exp(^/r)-l). 
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For the scalar field, M(^) in ([TJ is a matrix in the orbital momentum indices / 
and /' with matrix elements 

1 

"=\ 



M,j,{^)^d,{^R)J-^ £ K,„{2^L)H\,, . (5) 



Here, the function di{x) results from the T-matrix for the scattering on the sphere. 
For Dirichlet and Neumann boundary conditions on the sphere we note 

^rw=|^, ^rw = ^f^. (6) 

^vW [Kv{x)|^/x) 

In these formulas, R is the radius of the sphere, L is the separation between the plane 
and the center of the sphere, Iv{x) and A'v(x) are the modified Bessel functions. We 
introduced the notations v = Z + 1/2, v' = /' + 1/2 and v" = /" + 1/2, which will 
be used throughout the paper 

The factors H^^^ in (|5]l result from the translation formulas. Their explicit form is 

< = V(2/ + l)(2/' + l)(2/" + l)(^J^')(^j;^^'), (7) 

where the parentheses denote the 3y'-symbols. 

The above formulas are for Dirichlet boundary conditions on the plane. For Neu- 
mann boundary conditions on the plane we have to reverse the sign in the logarithm 
in dill or, equivalently, to change the sign of M. The trace in ([T]i, 

Tr= £ £, (8) 

is over the orbital momenta truncated at some /,„. Of course, the final expression 
appears for /,„ — > oo. 

For the electromagnetic field, the matrix M is in addition a matrix in the two 
polarizations. These correspond to the TE and the TM modes in spherical geometry 
and we can represent the corresponding matrix elements M/// as matrixes (2x2), 

with the factors 



_ ^■[l"il"+l)-l{l + l)-l'{l' + l)] - 2m^L 



which follow from the translation formulas for the vector field. The factors resulting 
from the scattering T-matrices are 



drix) 



lyjx) 
Kv{x) 



drix) 
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{Iv{x)^/x)' 



{K^{x)^)' 



(11) 



When inserting these expressions into ([T]) or (JUi, the trace must be taken also over 
the polarizations and the orbital momentum sum is restricted by I > max(l, \m\). 



3 The low temperature expansion 



Due to the Boltzmann factor in Q, the low temperature expansion emerges from 
the expansion, for (^ ^ 0, of 



M(^)=Mo + Mi(L^)'+M2(L^)2 + M3(L^)3 + .... 
The coefficients M,- = M, (p ) are dimensionless functions of the ratio 

R 



(12) 



(13) 



Inserting the expansion (fT2l i into the trace of the logarithm and keeping only the first 
two odd orders we get 



Trlnil-M{^)) =Ni{pn+N,{p){L^f 



with 



= -Tr 


(1-Mo)"'Mi 


: 






= -Tr 


(1-Mo)"'M3 






-Tr 


(1-Mo)"'Mi 


(1-Mo)"'M2 


-1- 


((1-Mo) ^M^y 



(14) 



(15) 



which are functions of p like the M,'s. 

It must be mentioned that inserting ( fTSl ) into dUi we interchange the orders of the 
Hmits r — > and /„, — s- °a. Below we will see in which cases this is justified and in 
which it is not. With the expansion (fTSl l, the low-T contributions to the free energy 
dill are 



-N,ip)LT^ + ^N3ip)L'T' 

The corresponding contributions to the force are 

Xf^ d ;r d(LM(p)) ^2 ^^ d(L^JV3(p)) ^4 

^ dL 6 dL 15 dL 



(16) 



(17) 



As it will turn out there is only one contribution ^ T^ to the force (Section 4.2, 
below) and in all other examples considered in this paper the low-T expansion starts 
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from T^. This is in agreement with the findings of QWGIOI . In order to compare the 
resuhs we expand ( fTTj i for small separation a = L — R, 



ATf={c2R' + CiaR^)r 



(18) 



where the coefficients C2 and ct, were introduced in the same way as in OWGIOI 



4 Results for hard boundary conditions on the sphere 

In this section we consider hard boundary conditions on both, the plane and the 
sphere. We start with the scalar field. Here we have Dirichlet (D) or Neumann (N) 
boundary conditions and we denote the four combinations by (X, Y), where X stands 
for the sphere and Y stands for the plane. For example, (D,N) denotes Dirichlet 
boundary condition on the sphere and Neumann boundary conditions on the plane. 
We remind that Neumann boundary conditions on the plane appear from reversing 
the sign of M. 



4.1 The case DD 



In this case we have a non-zero 



Ni=p., 



(19) 



which is independent on the truncation Im. This is the only case where one of the 
function A^i or A^3 considered in this paper does not depend on the truncation. 

The function A^i, ( fT9] l, delivers the T^-contribution to the free energy which was 
found in IIBP09I . It does not contribute to the force since the dependence on L drops 
out. The next-order contribution is N^ . It is a rational function of p . The orders of the 




Fig. 1 The functions Ni(p) for the case (D,D) (left panel) and N\{p) for the case (D,N) (right 
panel) for several values of the truncation /„, . The limit of small separation corresponds to p = 1 . 
The dashed line coiTesponds to /„, = 0, i.e., to the pure s-wave contribution. Already for /„, > 1 
there is nearly no dependence on /„,. 
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polynomials in numerator and denominator grow with the order /„, of the truncation. 
This is a general feature and holds for all functions A^i and A^3 considered below 
except for those which vanish. 

We display Nj,{p) as a function of p in fig.[T](left) for several /„,. It is seen that 
already for Z„, > 1 the curves cease to change. In this way the free energy and the 
force have a well defined limit for /„, -^ °°. The coefficients C2 and C3 defined in ( fTSl l 
are shown in Table l4.l] and it is seen that these fit well to those found in IIWGIOII . 

At large separation, i.e., for a small sphere, only the lower orbital momenta are 
on work. The function A^i is given by eq. iT% . the function A^3 by 

A^3 = ^P + ^p'-^p' + 0(p4). (20) 

In this way, the leading order temperature correction to the free energy and to the 
force do not depend on separation, while the subleading, ^ T'^, correction does. 



/„ 



■2.756 -3.748 -3.770 -3.772 -3.772 -3.772 -3.772 -3.772 -3.772 
■5.512 -2.910 -2.500 -2.429 -2.426 -2.427 -2.426 -2.425 -2.425 



Table 1 The values of the coefficients ct and C3 defined in eq. )18t for the case (D,D) for several 
values of the truncation /„,. The corresponding values found in IWGIOI . eq. (25), ai'e C2 = —3.96 
and C3 = —2.7. 



4.2 The case DN 

In this case, as in the previous one, the dominating contribution is A^i . However, now 
it depends on the order /,„ of truncation. The first two orders are 

^'(^)|'™^o - 2 + p ' ^'(^)l'-' - i6 + 8p+4p3 + p4 ■ (21) 

In fig. [T|(right panel), A^i is shown as function of p for several values of Z„,. It is seen 
that there is a rapid convergence for large Z„,. The sign is reversed as compared to 
the case (D,D) like for parallel plates. 

Because of the more involved dependence on p as compared to ( fT9l ). Ni con- 
tributes to the force. Hence, for these boundary conditions, according to ( fTTI ). we 
have a contribution ^ T^ to the force. 
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Fig. 2 The functions Af;(p) for the cases (N,D) (left panel) and (N,N) (right panel) for several 
values of the truncation /,„. The limit of small separation corresponds to p = 1. The dashed line 
corresponds to /„, = 0, i.e., to the pure s-wave contribution. Already for /„, > 1 there is nearly no 
dependence on /,„ . 



4.3 The cases ND and NN 

In these two cases, which have Neumann boundary conditions on the sphere, the 
contribution of A^i is zero for all /,„ and the expansion starts with A^3 . These functions 
share the common features discussed above. We displayed both cases in fig.|2] 

We mention that in IIBP091 we considered only the contribution from /„, == and 
missed the higher order terms. 

The expansions for large separation is 



A^- 



ND) 



A^- 



3 

(NN) 



24' 



192' 



24' 



384' 



+ 0(p"), 
0(p"). 



(22) 



The first two terms are the same (up to the sign), higher orders are different. 



4.4 The electromagnetic field with conductor boundary conditions 



For the electromagnetic field, the matrix elements M/// are given by eq. (|9]l. We 
expand them in powers of £, as before and obtain an expansion in parallel to eq. 



m. 



(^) 



(L^)'+M2(L^)2 + M3(L^)^- 



(23) 



It turns out that the matrixes M, are diagonal in the polarizations for / = and / = 2 
and anti-diagonal for i=l. Therefore in eq. ( fT4] i. the first contribution, A^i , vanishes 
and from Nj, only the term in eq. (fTST i does not vanish. From this structure it follows 
also that the contributions from the two polarizations do not mix in A^3 and we can 
consider these separately. Now the further calculations go in the same way as for 
the scalar field and we have calculated the functions Ni{p) for both polarizations. 
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We displayed them in fig. [3] While the TE mode gives a result similar to the scalar 




Fig. 3 The functions A'3 (p ) for the electromagnetic field with conductor boundary conditions for 
the TE polarization (left panel) and for the TM polarization (right panel) for several values of the 
truncation /„, . The limit of small separation corresponds to p = 1 . The dashed line corresponds to 
/m = 1, i.e., to the pure p-wave contribution which is the lowest one in the electromagnetic case. In 
the TE case, already for /„, > 1 , there is nearly no dependence on /„, . For the TM case, at p < 1 , 
there is no convergence for growing /,„. We displayed until /„, = 10. 



case in the sense that the limit /„, -^ 00 is approached very fast, the picture for the 
TM mode is different. Here we observe, for p close to unity, p < 1, contributions 
growing with /„, . This must be interpreted as a noncommutativity of the limits T ^ 
and l,n — > 00. As a consequence, at small separation, we have to expect contributions 
decreasing for T — j- slower than T"*. 

For large separation we find the following expansions. 



M 



TE 



1 



Nr 



3 
TM 



1 

12' 



2^16 



^pV0(p>'). 



According to (fTSl l. the leading order contribution to the free energy is 



337-4 



Tt'R^ 



At^ = ^R'T" -^'—T" 



15 



60 L3 



(24) 



(25) 



and we observe a T'^-contribution to the force (from the second term). The first term 
coincides with the corresponding term in eq. (6) in OCDNLRIOI while the second is 
beyond of what is displayed there. 



5 Results for a dielectric ball in front of a conducting plane 



For a dielectric ball the formulas of section 2 remain valid except for the functions 
dj^{x), eq. (fTTT i. These must be substituted by 
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JTE/„^ _ 2 ^/esi{x)s',{nx)-y^s'i{x)si{nx) 



driz) = - V V- ^ ^ (26) 

n \/zei [x)Si (nx) — y^ej {x)si [nx) 



with the refraction index n = ^/sjl and st{x) = \J %xl2ly{x) and eiix) = \flxj% 
K^ix) are the modified spherical Bessel functions. The function d™{z) can be ob- 
tained from ( [26l ) by interchanging e and jx. 

Inserting these formulas into (|9]l and calculating the entries in eq. (l23l ) we see that 
the structure of the matrixes M, remains the same. In this way we have a separation 
into TE and TM modes as before. As a consequence, we have only T^ contributions. 

Now we calculate the function A^3 (p ) for the case of a fixed £ and for an e taken 
from the plasma model. 



5.1 Fixed permittivity e 

For fixed e we consider two cases. First we put /i = 1. In this case it turns out that 
Nt,{p) =0 for the TE mode. This means that the corresponding low-T expansion 
starts from a higher power in T which we do not consider here. For the TM mode 
the function A^3 (p ) is shown in fig. |4] (left panel). It depends on the truncation. For 
e close to unity it stabilizes rapidly, for higher e slower 
At large distances, p ^- 0, we found 

In dilute approximation, e = 1+5, 5 ^ 1, only the lowest orbital momenta con- 
tribute until the order quadratic in 8, 

N]'' = ^p'8-^p\2 + p')5' + ff{5'), (28) 

higher orders are more complicated to obtain. 

As the second case we consider e = 1. Here the contribution of the TM mode to 
N2{p) is zero and we are left with the TE contribution. This function is very similar 
to that in the first case, however, different in details. It is shown in fig.|4](right panel). 
It stabilizes much faster when lifting the truncation as in the previous case. For large 
separations it reads 

^' -'W^f WTW^ "24(M + 2)3^ +^(^ ^- ^^^^ 

The difference (up to the sign) starts with order p^. In dilute approximation we 
found in the first two orders the same expression as in eq. (l28b with reversed sign. 
Differences show up starting from the third order 
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Fig. 4 The functions N™{p) for the dielectric ball (left panel) with /i = 1, e = 1.5 (lower curve), 
e = 2.3 and e = 10, and N]^{p) (right panel) and with e = 1, ^ = 1.8 (upper curve), ^ = 2.3 and 
^ = 10. The dashed line is for l„, = I. 



5.2 Plasma model permittivity 



The permittivity derived within the plasma model for metals is 



CO 



(30) 



where COp is the plasma frequency and e is taken for imaginary frequency ^ . In- 
serting ( l30b into ( l26b . the calculation runs as in the previous cases. The results are 
shown in fig.|5] 




Fig. 5 The functions N^ip) for the dielectric ball with dispersion <30b of the plasma model for 
Op = 1. For the TE mode (left panel) the convergence for /„, — )> oo is rapid for all separations. For 



the TM mode (right panel), there is no convergence for /,„ 
The dashed curve corresponds to /„, = 1 . 



> oo at small separation, i.e., for p < 1. 



For the TE mode, the curves cease to change already for /„, = 2 within the preci- 
sion of the plot. The analytic expressions are rational function of p and of hyperbolic 



functions of (Op. For small (Op we observe 



N^ 



(P) = 



45 



0{(ol). 



(31) 
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A different picture we observe for the TM mode. Here the truncation can be removed 
for p < 1 only. At close separation, i.e., for p < 1, the contributions do not tend to 
a finite limit for /,„ -^ °o. Again, we have to interpret this a non-commutativity of 
the limits f ^ and /„, -^ °°. Hence we expect for small separation lower powers in 
T. Also the behavior for cOp ^- is nonanalytic in the sense, that Ni{p)'^^ does not 
vanish in this limit, 

Nj^{p) = P^^±^+Oico„). (32) 

At large distances, the leading correction to the free energy for the TM mode 
does not depend on the plasma frequency, 

Nj^^^ + ^ip^). 
While for the TE mode the correction is sensitive to small plasma frequencies, 

TF /I 1 cothffflo) 
y3 (Op COp 

but saturates at 1/3 when (Op -^ °°. 



6 Conclusions 

In the forgoing sections we calculated the low temperature expansion of the free 
energy for a sphere or a dielectric ball in front of a plane. For the temperature 
dependent part At^ of the free energy we used the representation (|4]i involving 
the Boltzmann factor. Further we used a truncation of the orbital momentum sum, 
/ < /„„ and interchanged the limits T — > and /,„ — > °o. After that, the low-T expan- 
sion is obtained simply by expanding the matrices M into powers of ^ and taking 
the lowest odd one. In this way, the low-T expansion takes the generic form 

At.^^^2T^ + .^aT'^ + ---- (33) 

The coefficient ^2 is present for the cases (D,D), section 4. 1 (but independent on 
the separation) and (D,N), section 4.2. It is zero in all other cases where J^4 is the 
leading order contribution. 

For all examples considered in this paper, at finite separation, ^2 and ^4 have a 
finite limit for /,„ -^ °o. Hence the generic low-T behavior is given by eq. ( |33] |. This 
holds, for instance, at large separation. A different picture appears for small separa- 
tion, p ^- 1 . In some cases, the closer the separation, the worse the convergence for 
Im -^ °°- In these cases we do not have a result for T ^- 0. However, we can expect 
lower powers of T to appear This is a topic of future investigations. 
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